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Porous  composite  electrodes  play  a  critical  role  in  determining  the  performance  and  durability  of  solid 
oxide  fuel  cells,  which  are  now  emerging  as  a  high  efficiency,  low  emission  energy  conversion  technology 
for  a  wide  range  of  applications. 

In  this  paper  we  present  work  to  combine  experimental  electrochemical  and  microstructural  charac¬ 
terisation  with  electrochemical  simulation  to  characterise  a  porous  solid  oxide  fuel  cell  anode.  Using 
a  standard,  electrolyte  supported,  screen  printed  Ni-YSZ  anode,  electrochemical  impedance  spec¬ 
troscopy  has  been  conducted  in  a  symmetrical  cell  configuration.  The  electrode  microstructure  has  been 
characterised  using  FIB  tomography  and  the  resulting  microstructure  has  been  used  as  the  basis  for  elec¬ 
trochemical  simulation.  The  outputs  from  this  simulation  have  in  turn  been  compared  to  the  results  of 
the  electrochemical  experiments. 

A  sample  of  an  SOFC  anode  of  6.68  jim  x  5.04  [im  x  1.50  pim  in  size  was  imaged  in  three  dimensions 
using  FIB  tomography  and  the  total  triple  phase  boundary  density  was  found  to  be  1 3  purr  2 .  The  extracted 
length-specific  exchange  current  for  hydrogen  oxidation  (97%  H2,  3%  H20)  at  a  Ni-YSZ  anode  was  found 
to  be  0.94  x  1 0-10,  2.14  x  1 0~10,  and  1 2.2  x  1 0-10  A  purr1  at  800,  900  and  1 000  °C,  respectively,  consis¬ 
tent  with  equivalent  literature  data  for  length-specific  exchange  currents  for  hydrogen  at  geometrically 
defined  nickel  electrodes  on  YSZ  electrolytes. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Porous  composite  electrodes  play  a  critical  role  in  determining 
the  performance  and  durability  of  solid  oxide  fuel  cells  (SOFCs), 
which  are  now  emerging  as  a  high  efficiency,  low  emission  energy 
conversion  technology  for  a  wide  range  of  applications. 

Ni-YSZ  is  a  common  choice  for  the  anode  material  of  SOFCs  in 
both  electrolyte  supported  and  anode  supported  geometries  owing 
to  its  low-cost,  thermal  and  chemical  stability,  good  electronic  con¬ 
ductivity,  and  high  catalytic  activity  for  hydrogen  oxidation  [1]. 
Providing  intimate  contact  of  ionic,  electronic  and  pore  phases 
at  triple  phase  boundaries  (TPBs)  is  critical  to  optimising  electro¬ 
chemical  performance.  In  these  composite  porous  anodes,  TPBs  are 
distributed  throughout  the  electrode  bulk. 

Recently,  the  use  of  advanced  focused  ion  beam  (FIB)  tomo¬ 
graphic  techniques  has  enabled  direct  access  to  information  on  the 
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spatial  distribution  of  TPBs  in  an  SOFC  electrode  structure  (e.g., 
[2-5]).  The  so-called  TPB  density,  a  measure  of  the  abundance  of 
TPB  length  per  unit  volume,  is  an  important  metric  for  the  compar¬ 
ison  of  different  electrode  structures  as  it  is  generally  viewed  that  a 
greater  TPB  density  will  increase  electrode  performance  for  a  given 
materials  set  [6]. 

The  electronic  conductivity  of  Ni  is  typically  several  orders  of 
magnitude  higher  than  the  ionic  conductivity  of  YSZ  [7].  As  the 
electronic  conductivity  of  YSZ  is  negligible  under  SOFC  operating 
conditions,  the  transport  of  charged  species  will  occur  preferen¬ 
tially  in  the  Ni  phase  and  the  electrochemically  active  region  will 
be  located  in  a  region  close  to  the  electrode-electrolyte  inter¬ 
face;  Primdahl  and  Mogensen  [8]  estimate  the  active  thickness  for 
Ni-YSZ  anodes  to  be  under  20  pun.  Therefore  there  is  diminish¬ 
ing  benefit  in  continuing  to  increase  electrode  thickness  in  order 
to  maximise  the  number  of  triple  phase  boundaries  [9],  as  much  of 
the  added  thickness  will  be  electrochemically  inert.  However,  ques¬ 
tions  remain  regarding  the  optimum  thickness  for  a  given  electrode 
structure  and  material. 

In  composite  electrodes,  the  distribution  of  triple  phase  contacts 
throughout  the  electrode  bulk  allows  distributed  current  genera¬ 
tion  in  a  3D  volume.  As  the  fuel  cell  community  places  increasing 
emphasis  on  the  importance  of  microstructure  for  electrochemical 
optimisation,  it  is  essential  that  relevant  electrochemical  param- 
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Nomenclature 

fn  material  fraction 

I<n  conductivity  (S  m-1 ) 

0n  volume  fraction 

Vn  potential  (V) 

Ptpb  TPB  density  ( pun-2 ) 

R  gas  constant  (J  mol-1  K-1 ) 

T  temperature  (°C) 

F  Faraday’s  constant  (C  mol-1 ) 

j  charge  transfer  current  per  unit  length  (A  |jinrr 1 ) 

x  distance  from  interface  (boundary  condition) 

pH2  partial  pressure 

r]  over-potential  (V) 

P  transfer  coefficient  =  0.5 

Dh2  diffusivity  (m2  s-1 ) 

Fio  ionic  activation  energy  (J  mol-1 ) 

Ea  activation  energy  for  hydrogen  oxidation  (J  mol-1 ) 

y  pre-exponential  factor  (A  pan-1 ) 

cr j0  ionic  conductivity  (S  m-1 ) 

Rcy  charge  transfer  resistance  (A  cm-2 ) 

z  number  of  electrons  transferred 

EIS  electrochemical  impedance  spectroscopy 

TPB  three  phase  boundary 


eters  can  be  quantified.  A  key  measure  of  electrode  performance 
is  the  exchange  current  density  which  provides  a  measure  of  the 
propensity  of  a  given  electrode  material  to  catalyse  an  electrode 
reaction  at  a  given  temperature. 

In  order  to  effectively  compare  the  electrochemical  performance 
of  composite  electrodes,  exchange  current  density  should  be  con¬ 
sidered  per  unit  volume,  or  per  unit  length  of  TPB.  Obtaining  such 
data  is  difficult,  but  experiments  using  geometrically  defined,  dense 
Ni  structures  on  YSZ  substrates  [6,10]  have  provided  some  electro¬ 
chemical  data,  reported  per  unit  length  of  TPB. 

Modelling  of  fuel  cells  over  many  length  scales  from  fundamen¬ 
tal  electrochemistry  [11],  through  cell  level  [7,12,13],  to  system 
and  stack  level  [14-16]  provides  valuable  insight,  allowing  us  to 
explore  operation  under  steady  state  and  transient  conditions. 
Models  also  afford  easy  manipulation  of  variables  and  observa¬ 
tion  of  their  effects;  this  information  can  be  used  to  predict  and 
enhance  cell  performance.  Numerical  simulations  provide  valuable 
information  on  cell  stress,  potential,  current  density  and  temper¬ 
ature  in  electrodes  and  may  also  predict  the  effects  of  changes  in 
microstructural  parameters  on  electrode  percolation  and  perfor¬ 
mance  [7,17]. 

Models  that  explore  phenomena  at  the  cell  level  require 
information  on  the  microstructural  framework  so  that  the  electro¬ 
chemical  simulation  can  be  performed;  this  may  be  in  two  [18]  or 
three  [7,1 7]  dimensions.  The  generation  of  microstructures  that  are 
representative  of  real  systems  is  currently  the  subject  of  a  signifi¬ 
cant  research  effort. 

Numerical  generation  is  one  possible  route  to  obtaining  these 
frameworks,  for  example  using  Monte  Carlo  techniques  to  approx¬ 
imate  the  packing  of  the  three  constituent  phases  [19].  The 
complexity  and  accuracy  of  the  frameworks  derived  in  this  manner 
varies  greatly;  however,  mono-disperse,  homogenous  and  spher¬ 
ical  particles  are  common  assumptions  (e.g.,  [18]).  More  complex 
models  may  also  consider  the  effects  of  sintering  on  the  material 
microstructure  (e.g.,  [7,17]),  but  such  techniques  can  be  computa¬ 
tionally  expensive  and  are  constrained  by  a  limited  understanding 
of  the  evolution  of  actual  microstructures  during  sintering  pro¬ 
cesses. 


The  availability  of  accurate  information  concerning  the  3- 
dimensional  distribution  of  TPBs  in  an  electrode,  for  example  from 
FIB  tomography,  provides  an  opportunity  to  combine  experiment 
with  simulation.  By  mimicking  experimental  microstructures  com¬ 
putationally,  one  gains  the  ability  to  compare  experimental  results 
with  the  output  of  electrochemical  simulations,  thereby  offering 
a  potential  tool  for  model  validation.  In  turn,  this  will  improve 
the  accuracy  of  the  simulation  which  may  then  be  used  to  design 
optimum  electrode  microstructures. 

In  this  paper  we  present  work  to  combine  experimental 
electrochemical  and  microstructural  characterisation  with  elec¬ 
trochemical  simulation.  Using  a  standard,  electrolyte  supported, 
screen  printed  Ni-YSZ  anode,  electrochemical  impedance  spec¬ 
troscopy  has  been  conducted  in  a  symmetrical  cell  configuration. 
The  electrode  microstructure  has  been  characterised  using  FIB 
tomography  and  the  resulting  microstructure  has  been  used  as 
the  basis  for  electrochemical  simulation.  The  outputs  from  this 
simulation  have  in  turn  been  compared  to  the  results  of  the  elec¬ 
trochemical  experiments. 


2.  Part  1— experimental 

Ni-YSZ  symmetrical  cells  of  1 1  mm  diameter  have  been  pre¬ 
pared  by  screen  printing  commercially  available  NiO-YSZ  ink 
(Cerampaste  manufactured  by  Ceramtec)  onto  commercially  avail¬ 
able  pre-sintered,  tape  cast  YSZ  pellets  (Nextech  Materials).  The 
single  layer  screen  printed  anodes  were  then  sintered  at  1350  °C. 
Reduction  of  NiO  to  Ni  followed  the  manufacturers’  guidelines:  cells 
were  heated  to  600  °C  in  97%  N2,  3%  steam,  at  600  °C,  the  concen¬ 
tration  of  H2  was  increased  to  97%  over  a  period  of  1  h,  where  the 
cells  were  held  for  1  h  before  cooling  in  N2. 

Electrochemical  impedance  spectroscopy  (EIS)  has  been  used  to 
characterise  the  electrochemical  performance  of  the  cells  at  open 
circuit.  EIS  was  performed  in  a  symmetrical  cell  configuration  using 
an  Autolab  Impedance  Analyser  system  (PGSTAT30-Eco  Chemie). 
EIS  scans  from  lOOmFIz  to  1  MHz  were  carried  out  at  tempera¬ 
tures  ranging  from  800  to  1000°C  in  gas  environments  with  H2 
varying  from  48  to  97%  with  3%  steam,  at  a  total  gas  flow  rate 
of  100  mln  min-1.  Equivalent  circuit  analysis  was  performed  using 
ZView2  (Scribner). 

FIB  tomography  was  used  to  characterise  the  Ni-YSZ  electrode 
microstructure  after  electrochemical  characterisation  had  been 
performed.  A  Zeiss  xB-1540  FIB-FEG  SEM  was  used  in  this  study. 
The  ion  and  electron  beams  were  separated  by  52°  allowing  simul¬ 
taneous  milling  and  imaging  to  be  performed. 

Prior  to  imaging,  the  electrode  surface  was  carefully  prepared 
in  order  to  minimise  negative  effects  such  as  streaking  and  re¬ 
deposition.  By  maximising  the  quality  of  the  raw  images  obtained, 
it  has  been  possible  to  conduct  automated  image  analysis  of  the 
resulting  data  set.  To  do  this  the  entire  sample  was  sputter  coated 
with  gold  and  firmly  attached  to  a  microscopy  stub  using  silver 
paint,  thereby  providing  a  low  resistance  path  to  earth  to  dissi¬ 
pate  the  electronic  beam  current  and  prevent  charging.  Once  an 
area  of  interest  was  identified,  a  thin  layer  of  Pt  was  deposited  onto 
the  local  electrode  surface  using  a  standard  FIB  gas  injection  sys¬ 
tem;  the  Pt  layer  serves  to  minimise  streaking  or  curtaining  effects 
that  arise  due  to  the  incidence  of  the  ion  beam  on  a  non-uniform 
tomography  [20]. 

Using  high  beam  currents,  a  U-shaped  trench  was  then  milled 
around  the  deposited  Pt,  as  outlined  in  Fig.  1.  The  trench  pre¬ 
vents  milled  material  becoming  re-deposited  on  the  face  of  interest. 
Finally,  the  face  of  interest  is  polished  using  low-beam  currents  to 
remove  any  re-deposited  material. 

Through  optimisation  of  the  sample  preparation  procedure  for 
charge  dissipation,  and  by  using  improved  automated  image  anal- 
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Fig.  1.  Schematic  illustration  of  the  FIB  milling  technique,  along  with  a  typical  image  from  the  anode  section. 


ysis  procedures,  it  was  not  necessary  to  perform  the  lift-out  sample 
preparation  technique  described  in  our  previous  study  [4]. 

Using  the  Zeiss  xBl  540,  it  is  possible  to  collect  tomography  data 
in  two  ways,  firstly  by  milling  and  imaging  continuously,  capturing 
images  whilst  the  ion  beam  is  still  incident  in  the  area  of  interest. 
However  whilst  this  technique  has  the  advantage  of  speed,  it  is 
sometimes  possible  to  observe  signal  noise  due  to  the  interference 
of  the  ion  and  electron  beams.  The  second  method  is  to  periodically 
interrupt  milling  whilst  an  image  is  captured;  resuming  once  the 
electron  beam  scan  has  been  completed.  Both  processes  are  fully 
automated.  In  this  study  the  second  method  of  image  collection 
was  adopted. 

Milling  was  conducted  at  lOOpA  with  ion  beam  accelerating 
voltages  of  30  kV.  Low  voltage  imaging  was  conducted  below  2  kV 
to  maximise  the  contrast  differential  between  Ni  and  YSZ  [21  ].  The 
largest  e-beam  aperture  was  used  to  flood  the  electrode  surface 
with  electrons  providing  some  charge  equalisation  for  the  inci¬ 
dent,  positively  charged  ion  beam  [22].  The  low  accelerating  voltage 
also  minimises  the  electron  beam  penetration  volume  [23]  which 
minimises  the  out-of-plane  information  obtained  from  the  porous 
substrate. 

A  sequence  of  170  images  was  collected;  pixel  dimensions  of 
20  nm  have  been  recorded  with  an  inter-slice  thickness  of  15nm. 
From  this  sequence,  100  images  have  been  selected  and  segmented 
for  further  analysis.  The  images  have  been  collected  in  SE2  (sec¬ 
ondary  electron)  mode  and  demonstrate  a  good  contrast  gradient 
between  the  Ni  and  YSZ  phase;  using  the  SE2  mode,  twin  bound¬ 
aries  in  the  Ni  have  not  been  observed,  which  were  problematic  in 
other  studies  using  back  scatter  detectors  [4]. 

Image  analysis  has  been  performed  using  Matlab  and  Amira 
software.  Images  are  manually  aligned  with  reference  to  the  sam¬ 
ple  edges  (see  Fig.  1)  before  they  are  cropped.  An  edge  preserving 
smoothing  filter  is  applied  in  3D  before  the  images  are  thresholded 
based  on  the  image  histogram,  corresponding  to  the  secondary 
electron  yield.  This  segmentation  process  has  a  recurrent  artefact; 
material  along  the  outside  edge  of  Ni  grains  has  been  erroneously 
identified  as  YSZ.  This  problem  arises  at  the  edge  of  material  grains 
because  of  the  high  contrast  gradient  and  is  removed  using  mathe¬ 
matical  morphology  operations  in  Matlab  [24].  Further  information 
regarding  the  application  of  this  technique  and  its  accuracy  are 
provided  in  Ref.  [4]. 

Once  the  images  have  been  successfully  segmented,  they  are 
combined  into  a  single  3D  matrix  which  is  then  analysed  for  per¬ 
colating  three  phase  boundary  contact.  The  calculation  framework 
used  has  been  presented  by  Golbert  et  al.  [7],  and  the  application  to 
FIB  serial  section  data  is  also  presented  in  [4].  In  this  case  the  grid 
dimensions  used  are  dependent  on  the  voxel  volume  of  the  raw 
data. 


3.  Part  1— results  and  discussion 

Electrochemical  impedance  spectroscopy  was  carried  out  in 
symmetrical  cell  configuration  at  temperatures  from  800  to 
1000°C;  results  of  the  EIS  are  presented  in  Fig.  2. 

Fig.  2(a)  shows  the  EIS  response  for  varying  temperature  in  97% 
H2, 3%  steam,  on  a  molar  basis.  Fig.  2(b)  shows  the  effect  of  varying 
H2  concentration  from  48  to  97%  at  900  °C:  impedance  spectra  were 
collected  over  the  frequency  range  100  mHz  to  1  MHz. 

As  spectra  in  both  figures  suggest,  three  major  processes  can 
be  identified  at  high,  middle  and  low  frequencies.  The  lowest  error 
equivalent  circuit  fit  has  been  obtained  by  implementing  a  model 


Fig.  2.  (a)  EIS  scan  for  Ni-YSZ  anode  symmetrical  cell  at  800,  900  and  1000  °C  (97% 
H2, 3%  steam  diluted  with  N2)  and  (b)  EIS  scan  for  Ni-YSZ  anode  symmetrical  cell  at 
900  °C  in  varying  H2  concentration  (48-97%  H2). 
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Fig.  3.  Equivalent  circuit  analysis  for  electrode  processes. 


Table  1 

Area  specific  resistances  extracted  from  EIS  data. 


Temperature  [°C] 

h2(%) 

Khf  [^cm2] 

Rm  [£2 cm2] 

Klf  [^cm2] 

1000 

97 

0.08 

0.86 

0.16 

1000 

73 

0.09 

1.45 

0.21 

1000 

48 

0.09 

2.06 

0.27 

900 

97 

0.39 

1.37 

0.14 

900 

73 

0.39 

1.71 

0.20 

900 

48 

0.40 

2.14 

0.23 

800 

97 

0.80 

1.54 

0.17 

800 

73 

1.05 

2.19 

0.29 

800 

48 

1.06 

2.76 

0.29 

with  two  capacitors  and  one  constant  phase  element,  each  in  paral¬ 
lel  with  a  resistor,  as  indicated  by  Fig.  3.  Fitting  data  of  area  specific 
resistances  for  three  processes  are  collected  in  Table  1,  the  mean 
error  for  the  least  square  fitting  was  less  than  0.03%  per  arc. 

According  to  the  data,  the  high  frequency  process  is  greatly 
affected  by  temperature  and  seems  to  be  largely  independent 
of  hydrogen  concentration.  Such  behaviour  is  consistent  with  a 
thermally  activated  charge  transfer  reaction,  reported  to  be  a  two- 
electron  process  ([25,26]).  Two  other  processes  are  identified  at 
medium  and  low  frequency  and  are  shown  to  be  affected  by  hydro¬ 
gen  concentration,  and  less  sensitive  to  temperature,  especially  at 
the  lowest  frequency.  These  processes  are  associated  with  diffusion 
phenomena  [27,28]. 

From  the  FIB  tomography  sequence,  100  slices  have  been  seg¬ 
mented  and  analysed,  using  a  voxel  dimension  of  20  nm  x  20  nm  x 
15nm  corresponding  to  a  sample  size  of  6.68  pm  x  5.04  pan  x 
1.50  |jim.  Whilst  larger  reconstructions  have  been  conducted  (e.g., 
[3,4])  this  volume  is  representative  of  the  screen-printed  electrode 
which  has  a  thickness  of  5  p,m. 

Fig.  4  shows  the  reconstructed  anode,  which  has  a  total  sample 
size  of  6.68  pan  x  5.04  pm  x  1.50  pm. 


Within  the  co-ordinate  geometry  of  the  electrode,  a  contigu¬ 
ous  ionic  pathway  is  required  between  the  electrolyte  and  TPB, 
and  the  gas  and  electronic  phases  must  be  percolated  between  the 
TPB  and  current  collector/gas  bulk.  The  reconstructed  volume  has 
been  arbitrarily  labelled  in  the  3D  Cartesian  system  with  the  vol¬ 
ume  6.68  pm  x  5.04  pun  x  1.50  pun  corresponding  to  the  x,  y  and  z 
dimensions,  respectively.  The  reconstruction  encompasses  the  full 
electrode  thickness,  labelled  here  as  the  y-plane.  The  percolation 
of  the  constituent  phases  must  be  considered  in  the  y-dimension 
in  agreement  with  the  required  co-ordinate  geometry  discussed 
above.  However,  because  the  direction  of  milling  (see  Fig.  1)  is  not 
defined  relative  to  the  co-ordinate  geometry  of  the  electrode  in 
the  x-z  plane,  percolation  must  be  considered  in  both  the  x-y  and 
y-z  planes.  This  factor  can  be  used  to  assess  whether  the  sample 
is  sufficiently  large  by  considering  the  homogeneity  of  the  recon¬ 
structed  samples  relative  to  the  bulk  electrode;  whilst  the  total  TPB 
density  will  be  the  same  irrespective  of  the  direction  of  percola¬ 
tion,  a  homogenous  sample  should  demonstrate  similar  percolated 
TPB  density  in  both  the  x-y  and  y-z  planes.  Where  the  recon¬ 
structed  volume  is  too  small  relative  to  the  sintered  particle  size 
distribution,  percolation  will  be  restricted  in  one  plane  because  the 
enhanced  percolation  in  a  single  phase  will  prevent  percolation 
of  the  remaining  phases.  For  example,  a  large  Ni  grain  relative  to 
the  dimension  of  the  reconstructed  volume  will  promote  electronic 
conductivity,  but  by  virtue  of  its  high  volume  fraction  will  hinder 
gas  and  ionic  phase  percolation,  reducing  the  number  of  percolated 
TPBs. 

In  this  sample,  a  total  TPB  density  of  1 2.99  fxnrr2  was  calculated, 
of  which  8.15  p,m-2  is  percolated  in  the  x-y  plane  (percolation  of 
63%)  and  9.34  p,m-2  is  percolated  in  the  y-z  plane  (percolation  of 
72%).  This  broad  agreement  in  percolation  in  both  planes  lends 
weight  to  the  homogeneity  and  therefore  the  statistical  relevance 
of  the  reconstructed  volume. 


Fig.  4.  3D  reconstruction  of  Ni-YSZ  anode— red  corresponds  to  YSZ  phase,  blue  to  Ni  and  pore  space  is  void.  Total  sample  volume  is  6.680  p,m  x  5.040  |xm  x  1.50  fxm.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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In  Part  2  of  this  work,  the  application  of  the  experimentally 
obtained  structure  as  the  basis  for  an  electrochemical  simulation 
is  considered,  and  the  estimation  of  the  length-specific  exchange 
current  is  presented  (exchange  current  per  unit  TPB  length).  The 
length-specific  exchange  current  has  also  been  extrapolated  from 
literature  results  for  model  anode  geometries;  its  calculation  and 
comparison  are  also  presented  in  order  to  provide  a  validation  of 
the  experimentally  derived  structure. 

4.  Part  2— electrochemical  simulation 


The  boundary  conditions  are  as  follows:  it  is  assumed  that  the 
electrode/current  collector  interface  (ion  insulating)  is  at  x  =  0  and 
the  electrode/electrolyte  interface  (electron  insulating  and  imper¬ 
meable  to  gas)  is  at  x  =  L.  Ph2o,  and  Ph2,  oo  indicate  the  partial 
pressure  of  H20  and  H2  at  the  inlet.  All  the  boundaries  at  the  direc¬ 
tions  ofy  and  z  are  assumed  to  be  insulating. 

In  our  model,  charge  transfer  is  assumed  to  be  the  rate  deter¬ 
mining  electrode  process,  and  thus  the  Butler-Volmer  equation  is 
used  to  calculate  the  current  density  as  shown  by  Eq.  (6),  where 
P  =  0.5  and  jo  is  exchange  current  per  TPB  length  in  A  |xnrr 1 . 


Transport  phenomena  in  porous  electrodes  are  intimately 
linked  to  electrode  microstructure;  specifically  to  the  porosity, 
the  pore  size  and  the  tortuosity  factor.  In  order  to  capture  this 
behaviour,  the  transport  phenomena  are  modelled  for  each  phase 
separately  and  linked  at  the  percolated  TPB,  where  charge  trans¬ 
fer  occurs.  For  transport  modelling,  the  microstructure  is  divided 
into  voxels,  with  each  voxel  represented  as  a  cube  of  pixels.  Each 
voxel  contains  different  volumes  of  the  three  phases.  Using  the 
VOF  methodology  [29],  the  phase  fractions  in  a  voxel  i  are  rep¬ 
resented  by  j/)Pore,  jj-ti0,  and  fi  eh  for  the  pore,  ionic  and  electronic 
phases,  respectively.  The  characteristics  of  a  voxel  i,  for  example  its 
conductivity  /c/,  are  assumed  to  be  dependent  on  the  characteristic 
values  in  the  pure  phases  (e.g.,  /c°,  n  =  pore,  io,  el)  and  the  phase 
fractions  in  that  voxel: 

Kj  =  Y&un  (!) 

n 

Note  that  materials  with  mixed  conductivity  can  be  handled  in  this 
framework  (for  example,  anodes  using  CGO  rather  than  YSZ),  since 
each  material  can  have  a  different  conductivity  for  each  charged 
species.  In  this  paper,  we  assume  that  each  phase  conducts  only 
one  species.  Assuming  steady  state  and  ideal  gases  in  the  pores, 
and  considering  only  the  diffusion  of  H2  and  H20,  we  model  the 
ionic  and  electronic  potentials  in  an  SOFC  anode  using  the  following 
system  of  equations: 


^(^io/io^^io)  —  PTPR/ 

(2) 

^(ffei/eiVVei)  =  2prpBj 

(3) 

v  (D“$’°"VPh,o) 

=  PTPBi 

(4) 

v(-DH/„tv  \ 

~^pBi 

(5) 

where  ptpb  is  the  local  TPB  density  (in  [xm  [xnrr3 )  and  j  is  the  charge 
transfer  current  per  unit  TPB  length,  crio  and  crel  are  the  ionic  con¬ 
ductivity  and  the  electronic  conductivity,  respectively,  DH2o  and 
Dh2  are  the  diffusion  coefficient  of  H20  and  H2,  respectively, /io,/el 
and  /pore  are  the  volume  fraction  of  ion  conducting,  electron  con¬ 
ducting  and  pore  phases,  Vi0  and  Ve\  are  the  ionic  potential  and  the 
electronic  potential,  respectively,  Ph2o  and  Ph2  are  the  partial  pres¬ 
sure  of  H20  and  H2,  Pis  the  Faraday  constant,  R  is  the  gas  constant, 
and  T  is  the  temperature.  All  quantities  expect  for  F  and  R  depend 
on  the  position  (x,  y,  z). 

These  phenomena  are  modelled  in  each  voxel  i,  and  the  volume- 
averaged  values  are  taken  over  all  the  voxels.  When  integrating  Eqs. 
(2)-(5),  a  small-enough  step  of  1  x  10-8  is  taken  in  the  variables 
for  finite  differencing  so  as  to  avoid  errors.  Vio  can  be  integrated 
from  0  V  up  to  1 .23  V,  depending  on  the  operating  temperature.  The 
electronic  conductivity  of  Ni  (cxel)  is  assumed  to  be  high  enough  to 
neglect  its  electronic  resistance,  so  the  electrode  potential,  V^,  is 
assumed  to  be  zero  throughout  the  electrode.  Both  Ph2o  and  Ph2 
can  be  integrated  over  the  range  0-1. 


j=j  o 


-  exp 


2(1  -  P)F 
RT 


(6) 


For  the  SOFC  anode  there  is  a  relation  for  the  potential  terms  as 
below: 


Veq  -  r)  =  Vio  -  Vel  (7) 

where  rj  is  the  over-potential,  Vi0  the  ionic  potential,  and  Ve\  the 
electronic  potential,  and  Veq  is  the  potential  at  equilibrium,  defined 
by  Eq.  (8): 


RT 

Veq  =  \Z0e  —  2/ In 


(8) 


V0e  is  the  standard  potential  which  is  determined  as  a  function 
of  the  Gibbs  free  energy  change  involved  in  the  reaction,  taking 
the  partial  pressures  of  H2,  H20  and  02  to  be  at  standard  pressure. 
The  variation  of  V0e  over  the  temperature  range  of  800-1 000  °C  is 
calculated  using  the  following  equation  [30]: 

Voe  =  1.23  -  2.304  x  10“4(T  -  298.15)  (9) 


Combining  Eqs.  (7)  and  (8)  gives  an  expression  for  the  over¬ 
potential: 


ri  =  Ve q  -  (V,0  -  Vel)  =  Voe  +  Vel  -  Vio  -  —In 


RT,  [Ph2o 


Ph, 


RT 

=  Voe  +  Vel-Vio  +  —\n 


(10) 


Substituting  (10)  into  (6)  results  in  a  definition  of  the  charge 
transfer  current  as  a  function  of  the  state  variables: 


r  f  ,  j 

j  =  Jo 

exp 

W (Voe  +  Ve,  -  Vfo) 

Ph2 

Ph2o 


-exp 


^(  Voe  +  Vel  -  Vi0 


Ph2o 

Ph2 


(11) 


The  ionic  conductivity,  (Tio,  is  estimated  by  using  Eq.  (12)  [31  ]: 


where  cr?  =  3.6  x  107  S  nrr1  is  the  pre-exponential  factor  (or  refer¬ 
ence  ionic  conductivity)  and  Fio  =  8.0  x  1 04  J  mol-1  is  the  activation 
energy  for  the  ionic  transport  [32]. 

The  system  of  equations  represented  by  Eqs.  (2)-(5),  (10),  (11) 
and  the  boundary  conditions  in  Table  2  is  solved  by  developing 
mass/charge  balances  for  each  voxel  surrounding  the  discretization 
points,  using  the  finite  volume  method  [29].  The  details  of  the  solver 
methodology  can  be  found  elsewhere  [7].  The  input  parameters 
used  are  given  in  Table  3. 

The  exchange  current  per  unit  TPB  length,  j,  is  fitted  using  an 
iterative  approach.  The  current  density-over-potential  response 
predicted  by  the  electrochemical  model  for  the  given  SOFC  anode 
microstructure  and  a  given  value  of  the  exchange  current  per  unit 
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Table  2 

Electrochemical  simulation  boundary  conditions. 
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@x  =  0 

@x  =  L 

dVioldx  =  0 

14,  =  0—1.23 

Vd  =  0 

dVel/dx  =  0 

Ph20  —  Ph20,  OO 

o 

II 

£ 

o 

if 

Ph2  =  Ph2  ,  oo 

dPH2/dx  =  0 

Table  3 

Electrochemical  simulation  input  parameters.  In  this  work,  length-specific  exchange 
current  is  used  as  a  fitting  parameter  (see  text). 


Temperature  (°C) 

800 

900 

1000 

Ionic  conductivity  (S  |xm-1 ) 

4.28  x  10“6 

8.40  x  10“6 

14.7  xl0“6 

Electronic  conductivity 
(Spjrr1)  [33] 

2.4 

2.17 

1.85 

Dh2  (m2  s”1)  [34] 

4.11  xl0“5 

4.11  x  10“5 

4.11  x  10“5 

DH2o(m2  s-1)  [34] 

4.11  x  10“5 

4.11  x  10“5 

4.11  x  10“5 

Voe  (V) 

1.05 

1.028 

1.005 

x  1CT 


1.2- 


c 
1 
3 

&  0.8b 


IJJ  0.6 


W  0.4 


0.2 


0L 

950 


+ 

model 

+ 

♦ 

Electrode  B 

• 

Electrode  C 

X 

Electrode  D 

- 

♦ 

1000 


1050 


1100  1150 

Temperature  (K) 


1200 


1250 


1300 


length  is  used  to  obtain  a  predicted  area  specific  resistance  for 
charge  transfer.  This  value  is  compared  to  the  area  specific  resis¬ 
tance  for  charge  transfer  as  estimated  from  EIS  measurements  (ftH f 
in  Table  1).  The  exchange  current  per  unit  length  is  adjusted  until 
good  agreement  is  obtained. 

The  following  procedure  is  used  to  compute  the  area  specific 
resistance  for  charge  transfer.  At  low  over-potentials,  it  is  possible 
to  assume  straight  line  behaviour  for  the  current-voltage  response 
according  to  the  low-field  Tafel  approximation  [8],  which  assumes 
charge  transfer  limited  (Butler-Volmer)  kinetics.  The  gradient  of 
the  voltage-current  density  plot  at  low  over-potential  (rj  <  20  mV) 
therefore  defines  the  charge  transfer  resistance. 

By  performing  the  simulation  procedure  outlined  here,  the 
exchange  current  per  unit  TPB  length  for  a  Ni-YSZ  system  in 
97%  H2,  3%  steam  has  been  calculated  at  three  different  temper¬ 
atures.  The  length-specific  exchange  current  is  estimated  to  be 
9.4  x  10"11,  2.14  x  10-10  and  1.22  x  lO^Afxm-1  at  800,  900  and 
1000°C,  respectively.  Fig.  5  shows  the  results  of  this  iterative  pro¬ 
cedure  at  three  temperatures  at  which  the  area  specific  resistances 
are  0.803  £2  cm2,  0.387  £2  cm2  and  0.076  £2  cm2,  respectively,  cor¬ 
responding  well  to  the  experimentally  measured  charge  transfer 
resistances  shown  in  Table  1,  i.e.  ]?hf-  This  part  of  the  EIS  response 
is  appropriate  for  comparison  as  our  electrochemical  model  cur¬ 
rently  considers  only  electron  transfer;  we  do  not  yet  account  for 


Current  density  (A/cm2) 


Fig.  5.  Over-potential  as  a  function  of  current  density  (normalised  for  geometri¬ 
cal  electrode  area)  as  predicted  by  electrochemical  simulation  (data  points  from 
electrochemical  simulation,  straight  line  approximation  provides  ASR). 


Fig.  6.  Calculated  variation  of  length-specific  exchange  current  with  temperature 
for  different  geometrically  defined  electrodes  at  open  circuit. 


additional  contributing  processes  such  as  adsorption  and  surface 
diffusion,  which  will  be  considered  in  future  studies. 

In  principle,  the  fitted  values  of  the  exchange  current  density  can 
be  compared  to  experimental  data  obtained  using  dense,  geometri¬ 
cally  defined  Ni  electrodes,  deposited  on  a  smooth  YSZ  surface.  This 
enables  the  isolation  of  triple  phase  contact  points,  the  length  of 
which  may  be  subsequently  measured.  In  doing  so,  electrochemical 
measurements  can  be  directly  correlated  with  TPB  length.  However 
there  are  very  few  such  data  sets  reported  in  the  literature.  Bieberle 
et  al.  [6]  have  deposited  dense  Ni  by  magnetron  sputtering  onto  a 
YSZ  single  crystal  surface.  The  deposited  Ni  was  subsequently  pat¬ 
terned  using  photo-lithography  to  alter  the  effective  TPB  length. 
Four  model  electrode  samples  were  prepared,  with  TPB  lengths 
ranging  from  0.04  to  8.738  m  (labelled  here  as  electrodes  A-D). 
Using  data  reported  by  Bieberle  et  al.,  it  has  been  possible  to  extrap¬ 
olate  data  on  the  length-specific  exchange  current;  the  calculation 
method  is  outlined  here. 

Fig.  1 1  from  Bieberle  et  al.  [6]  shows  the  variation  of  polarisation 
resistance  with  TPB  length.  For  the  geometrically  defined  system, 
it  is  assumed  there  is  no  diffusion  limitation  and  the  resistance  is 
therefore  attributed  to  charge  transfer  only.  Using  the  following 
equation  it  is  possible  to  extract  the  exchange  current  density  from 
the  polarisation  resistance  assuming  Butler-Volmer  kinetics.  This 
allows  direct  comparison  with  the  length-specific  exchange  current 
density  obtained  from  the  model  assuming  Butler-Volmer  kinetics: 


RT 

~  z~-  Rcr  ■  F 


(13) 


The  temperature  dependence  of  exchange  current  per  unit 
length  is  extrapolated  using: 


Jo  =  Y  ■  exp 


(14) 


where  EA  is  the  activation  energy  and  z  is  the  number  of  elec¬ 
trons  transferred.  The  pre-exponential  factor  is  calculated  from 
the  experimental  data  obtained  at  800  °C.  Literature  values  for  the 
activation  energy  of  hydrogen  oxidation  at  nickel  electrodes  [32] 
were  used  to  extrapolate  the  behaviour  of  exchange  current  with 
temperature,  this  is  shown  in  Fig.  6.  By  extrapolating  the  changes 
in  length-specific  exchange  current,  it  is  possible  to  compare  the 
experimental  results  with  our  predictions;  this  is  shown  in  Fig.  6 
and  Table  4. 

In  Table  4,  the  exchange  current  per  unit  TPB  length  extrapo¬ 
lated  from  measured  data  (electrodes  A-D)  is  presented  alongside 
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Table  4 

Temperature  profile  for  length-specific  exchange  current:  experimental  data  from 
Bieberle  et  al.  (electrodes  A-D)  and  predictions  from  our  model. 


Electrode 

Jo  (A  fxm-1 

)  at  temperature  (°C) 

800 

900 

1000 

A 

4.69E-10 

1.22E-09 

2.73E-09 

B 

1.07E-10 

2.77E-10 

6.21E-10 

C 

1.20E-10 

3.13E-10 

7.01E-10 

D 

1.21E-10 

3.14E-10 

7.02E-10 

Model 

0.94E-10 

2.14E-10 

12.2E-10 

exchange  current  data  predicted  using  the  simulation  presented 
here  (labelled  “model”).  The  predicted  data  can,  in  most  cases,  be 
seen  to  compare  favourably  with  length-specific  exchange  current 
data  extrapolated  from  the  literature.  The  main  difference  that 
can  be  observed  is  for  electrode  A,  which  may  reflect  the  small 
TPB  length  of  this  experimental  sample  (0.04  m)  compared  with 
the  other  patterned  electrode  samples  (ranging  between  2.5  and 
8.7  m),  which  would  be  expected  to  increase  experimental  error. 
Consequently  the  value  for  electrode  A  was  not  used  further  in 
our  analysis,  but  an  average  of  electrodes  B-D  was  used,  namely 
1.16  x  10-10,3.01  x  10-10, and 6.75  x  lO”10  Apum-1  at800,900and 
1000  °C,  respectively. 

Furthermore  we  note  that  Utz  et  al.  [10]  have  recently  sug¬ 
gested  that  for  patterned  Ni  electrode  geometries  on  YSZ  substrates, 
increases  in  TPB  length  of  up  to  30%  can  be  observed  when  heating 
the  sample  from  room  temperature  to  SOFC  operating  temper¬ 
atures.  This  requires  further  study  to  refine  these  challenging 
measurements. 

It  is  also  evident  from  Table  4  that  agreement  between  our 
extracted  data  and  the  experimental  data  is  better  at  lower  tem¬ 
perature.  This  is  perhaps  not  surprising  as  we  have  assumed  charge 
transfer  to  be  rate  determining,  which  will  be  more  likely  at  lower 
temperatures,  as  evidenced  by  the  higher  charge  transfer  resistance 
reported  for  lower  temperatures  in  Table  1. 

5.  Conclusions 

The  potential  to  link  electrochemical  and  microstructural  char¬ 
acterisation  with  electrochemical  simulation  has  been  explored. 
The  estimation  of  length-specific  exchange  currents  has  been 
demonstrated  using  a  combined  modelling  and  experimental 
approach,  which  has  compared  favourably  with  the  limited  exper¬ 
imental  data  available  from  the  literature. 

A  sample  of  an  SOFC  anode  of  6.68  p,m  x  5.04  [xm  x  1.50  [xm 
in  size  was  imaged  in  three  dimensions  using  focused  ion  beam 
techniques  and  the  TPB  density  was  found  to  be  12.99  p,m-2. 
This  same  anode  was  measured  using  electrochemical  impedance 
spectroscopy,  and  the  charge  transfer  resistance  determined  at 
temperatures  of  800,  900  and  1000°C. 

An  electrochemical  anode  model  has  also  been  used  to  represent 
the  system  using  the  imaged  anode  microstructure.  The  perfor¬ 
mance  predicted  by  the  model  has  been  compared  to  the  measured 
electrochemical  performance,  and  the  length-specific  exchange 
current  varied  until  agreement  was  reached  between  the  predicted 
and  measured  performance,  with  the  assumption  that  charge  trans¬ 
fer  was  rate  determining.  The  extracted  length-specific  exchange 
currents  for  hydrogen  oxidation  (97%  H2,  3%  FI20)  at  a  Ni-YSZ 
anode  were  0.94  xlO"10,  2.14  xlO-10,  and  12.2  x  1CH0  A  parr1 
at  800,  900  and  1000°C,  respectively.  These  values  are  consis¬ 
tent  with  equivalent  literature  data  for  length-specific  exchange 
currents  for  hydrogen  oxidation  at  geometrically  defined  nickel 
electrodes  on  YSZ  electrolytes;  namely  1.16  x  10-10,  3.01  x  10-10, 
and  6.75  x  1 0-10  A  |jnrr 1  at  800, 900  and  1 000  °C,  respectively,  not¬ 


ing  that  the  increasing  discrepancy  between  these  two  data  sets 
with  temperature  may  reflect  the  simplifying  assumption  made  in 
this  work  that  charge  transfer  is  rate  determining. 

Future  work  is  focused  on  fully  validating  the  electrochemi¬ 
cal  simulations,  with  the  aim  of  developing  a  powerful  tool  for 
electrode  design.  Furthermore,  it  should  become  possible  to  com¬ 
pare  the  results  of  3D  reconstructions  of  real  electrodes  with 
Monte  Carlo  simulations  of  electrode  microstructures,  allowing 
us  to  improve  the  accuracy  of  numerically  generated  electrode 
analogues.  Work  is  also  underway  to  expand  the  existing  elec¬ 
trochemical  simulations  to  extend  the  electrode  kinetic  model  to 
include  adsorption  and  surface  diffusion  processes. 
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